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We present a calculation of neutrino pair bremsstrahlung due to electron-electron scattering in a 
relativistic degenerate plasma of electrons. Proper treatment of the in-medium photon propagator, 
i.e., inclusion of Debye screening of the longitudinal part and Landau damping of the transverse part, 
leads to a neutrino emissivity which is several orders of magnitude larger than when Debye screening 
is imposed for the tranverse part. Our results show that this in-medium process can compete with 
(^-^ , other sources of neutrino radiation and can, in some cases, even be the dominant neutrino emission 

• mechanism. We also discuss the natural extension to quark-quark bremsstrahlung in gapped and 
' ungapped quark matter. 

(N : 

^ , PACS: 26.60.+C, 95.30.Cq, 97.60.Jd 

<: 

I I. INTRODUCTION 

(N . 

• For many astrophysical objects (e.g., supernovae or neutron stars) in which high densities and/or temperatures are 
^ ' attained, energy losses by neutrino emission are the driving force in their evolution [1-3]. The simplest process of /3- 

. and inverse /3-decay of nucleons, the so-called direct Urea process, leads to prodigious energy losses [4] but is possible 
' only in cases where free nucleons are present and, even in this case, can be kinematically forbidden. Similar direct 
Urea processes are also possible when hyperons [5] or quarks [6] are present and the associated emissivities dominate 
the cooling, but the existence of these exotic phases is subject to large uncertainties. In many cases, neutrinos can 
only be emitted through higher order processes in which many-body effects in the medium become important. A 
classic example is the plasmon decay process of neutrino pair emission where the very existence of plasmons is due 
to the medium's response to electromagnetic perturbations. Neutrino radiation via bremsstrahlung comprises a large 
'q^. subset of all neutrino emission processes where scattering of particles, through electromagnetic or strong interactions, 
I ' allows transfer of energy-momentum accompanied by the emission of a Z*^ boson which immediately decays into a v-V 
pair. Typical and well studied examples are, among others, the case of nucleon-nucleon, electron-ion and electron- 
^ phonon scattering in the neutron star's core and crust (see, e.g., [7] for a recent review). However, neutrino pair 
bremsstrahlung emission in collisions of relativistic particles has received little attention and only a few attempts 
exist in the literature, e.g., [8]. In this paper, we will be concerned with the importance of the medium in the 
modulation of the scattering interaction of relativistic electrons through polarization effects, and its impact on the 
neutrino emissivity. 
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In a QED plasma, static electric fields are screened while magnetic fields suffer damping only at non-zero frequency. 
For scattering in the relativistic limit, when magnetic effects are important, soft photon exchange can lead to a 
divergence from the unscreened magnetic interaction. To avoid this divergence, in early works on transport coefficients 
in a QED plasma, Debye screening was imposed in both polarization channels of the photon propagator in the covariant 
gauge, even though it only apphes to its longitudinal "electric" part. It was later duly noted [9,10] that Landau 
damping of the transverse "magnetic" part can generate a finite result in these calculations. Landau damping instead 
of Debye screening of the magnetic interaction implies higher scattering rates and hence smaller thermal and electrical 
conductivities and shear viscosity [10]. Similarly, higher scattering rates should result in larger neutrino emissivities 
for bremsstrahlung processes, and this is the issue we address here by calculating the neutrino pair emission from 
electron-electron bremsstrahlung in a degenerate relativistic electron plasma. Our results confirm these heuristic 
arguments and show that erroneously assuming Debye screening of the magnetic interaction can underestimate the 
neutrino emissivity by several orders of magnitude, particularly at the low temperatures of relevance to a variety of 
astrophysical environments. 
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This paper is organized as follows: in section II, we discuss the many-body effects that are relevant for the temper- 
ature and density region under consideration and which must be taken into account. Section III outlines the principal 
steps in the calculation of the neutrino bremsstrahlung emissivity from electron-electron collisions. We present a 
comparison of the results with other neutrino rates in three different contexts : from gapped and ungapped quark 
matter in strange stars in section IV, from hot matter in a black-hole accretion torus in section V, and finally in the 
more conventional medium of a neutron star, both in its crust and in its core, in section VI. Our conclusions and 
outlook are presented in section VII. Some numerical details are explained in the appendix. 

II. NEUTRINO BREMSSTRAHLUNG IN ELECTRON-ELECTRON COLLISIONS 

We are interested primarily in the effects of Landau damping on the bremsstrahlung rate since the screened case 
has already been considered in [8]. A unified treatment, with consistent inclusion of plasma effects, of neutrino 
bremsstrahlung in electron-ion and electron-phonon collisions and plasmon decay in the crust of neutron stars is 
described in [11]. Here, we derive the emissivity for bremsstrahlung in scattering of relativistic electrons, which is 
relevant to a host of astrophysical conditions, and the results are also extended to quark-quark bremsstrahlung, where 
the exchanged gluon is similarly dressed. Neutrino bremsstrahlung is computed from the set of tree-level Feynman 
graphs displayed in Fig. 1. Before performing the explicit computation, we discuss some many-body effects that can 
potentially be important for our temperature and (electron or quark) density regime. 

a) Landau damping: In the electron plasma, the exchanged photon is dressed by interaction with the electrons. 
The longitudinal photon is Debye screened, while the transverse photon is Landau damped. At the one-loop level, an 
all-order resummation leads to a photon propagator that is modified from its "bare" expression as follows: 

where K'^ = ui^ — and P^'^ , are the longitudinal and transverse projection operators in [12]. The longitudinal 
and transverse self-energies 11^ and Wt are in general complex functions of lo and k. These self-energies can be 
obtained in the hard dense loop approximation ^. Since the exchanged photon also has small energy, we will use 
approximate forms [12] in the nearly static limit w/k <^ 1: 

Hl = mjj and Ht = m\[ = i ^mf, . (2) 

For degenerate electrons, characterized by a chemical potential /Xe, 

m]} = e^lJ-l/^ , (3) 

but for the sake of generality, we use m,£i to denote the phase-dependent screening mass. Its explicit expression is to 
be understood by the particular case mider consideration. 

b) Vertex renormalization: Neutrinos are emitted with average energies that are thermal, and in the degenerate 
regime T/^^, <^ 1. This implies that their wavenumber k^, is such that k'^/mj-, <^ 1. Physically, this means that 
at sufficiently low temperatures, and large electron densities, the neutrino's wavelength is comparable to or exceeds 
the Debye charge radius. In this case, as shown by Leinson [11], the effective weak charge seen by the neutrino is 
renormalized. The correction factor due to this effect can be absorbed simply into "effective" electroweak couplings, 
and we take this into account in our calculations. 

c) Non-Fermi liquid effects: In normal quark matter with Nf quark flavors at a common chemical potential n, the 
absence of screening for the chromomagnetic gluon leads to a breakdown of Fermi liquid behaviour for quarks at a 
scale T ~ m exp(— 97r/4Q!c) [13] where rn^ = Nf-ac^'^ /tt is the electric screening mass, with ac = .9^/47r where g is the 
QCD coupling constant. The neutrino emission rate for the quark direct urea process is consequently enhanced [14]. 
This effect applies also to the electron plasma, but we ignore it here since the smallness of the QED coupling constant 
implies that this effect sets in only at exponentially small temperatures. 



^We are ignoring the smaller contribution from finite temperature efi'ects (hard thermal loops) since they are small for a 
degenerate QED plasma. 
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FIG. 1. The Feynman diagrams for neutrino bremsstrahlung in electron-electron scattering. The first 4 diagrams, enclosed 
within square brackets, are the "direct" contributions while the interchange (3 ^ 4) generates the "exchange" diagrams. The 
blob represents the one- loop resummed photon propagator (see text for details). 

d) LPM effect: Suppression of soft radiation due to multiple scatterings between electrons (the Landau- 
Pomcranchuk-Migdal or LPM effect) becomes important when the formation time t j of the neutrino wave-packet 
is comparable to the mean collision time tc between electrons. The ratio tf/tc, for a degenerate plasma (T <C /Ue) 
behaves approximately as sfoiTj /Zg if only screening is considered, and as a when damping is taken into account {a 
is the fine structure constant). In cither case, the ratio is much smaller than 1, so that an inclusion of the LPM effect 
would make only a 1% correction to our results. The large degeneracy at the Fermi surface implies that the mean 
time between electron-electron collisions is much greater than the formation time of the neutrino wave-packet. 



III. NEUTRINO EMISSIVITY FROM THE BREMSSTRAHLUNG PROCESS 



The neutrino emissivity from the bremsstrahlung process is 
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The subscripts i = 1 to 4 on p refer to electrons (1,2 label the incoming states and 3,4 the outgoing states), while i,f 
refer to total initial and final three- momentum/energy accordingly. Qi = {w\,q\),Q2 = (w2,?2) are respectively the 
neutrino and anti-neutrino four-momenta with Q = {Qa, q) = {lui +(-^2, qi+q2)- The symmetry factor s = 2. The phase 
space for electrons is convolved with the appropriate Fermi distribution functions npiEi) = l/(e(^*^'''=)/^ -|- 1) (in 
fcs = 1 units) and np = l — np, respectively. The matrix element M is obtained by the application of Feynman rules 
to the diagrams shown in Fig. 1. Since we work in the degenerate regime, the emissivity is dominated by contributions 
from low energy neutrinos (typically of order the temperature T, where T <C i^e) and Low's theorem [15] may be 
used to simplify the calculation of Qee-Br- As a consequence, the propagator of the internal fermion line can be 
approximated as 



Seip±Q) = ± 



Qo 



(5) 



A further approximation can be made by noting that interference between the direct and exchange diagrams is strongly 
suppressed relative to their squares taken individually, and hence may be neglected [8]. Then, in the low-energy limit. 
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the emission from the exchange diagrams simply doubles the overall result for the emissivity from the direct diagrams. 
Since there is a symmetry factor of s = 2 in the denominator of the emissivity expression eqn. (4), we may ignore the 
identity of the particles altogether and evaluate the emissivity from the direct diagrams alone (see caption Fig. 1). 

Aside from the equation for the dressed photon propagator (eqn.(l)), we also require the weak interaction Hamil- 
tonian (in the 4-Fermi form) given by 

= ■4'el\{cv - IbCAji'e ■ (6) 

The weak vector and axial- vector couplings are 

c\r = 1 + ^sin^Ow, c'X = l ■ 
4'^ = -l + 4sin2(?„,, c^'" = -l . (7) 

where the Weinberg angle sin^Ow = 0.231. Neutrino-pairs of all three flavors {e,ii,T) contribute incoherently to the 
emissivity, so that we may define a "total" squared couphng constant = ^^{{cy)^ + (cj^)^) = 6.696. The vertex 
renormalization effect discussed in the previous section introduces a suppression factor s = 0.449. This implies a total 
"effective" squared coupling constant c^g = 3.006. We will use the symbols cy and ca for the following calculation, 
and substitute the effective coupling constant only at the end. Evaluating the spin-summed squared matrix element 
with the low-energy approximations for neutrinos, and utilizing symmetries between the electron "labels" (1,2,3,4), 
we obtain (neglecting terms of 0{m1/ nD) 

IMP 2^fip^r^^ 

Wo 



s 

{Tl).ySnp = [{5<5p(Pi/xP37 + 1 ^ 3) g^),{pipP3s + 1^3)+ Pi- Psig^pg^s)} - (7 ^ <5)] , (9) 
L{Qi,Q2r^ = miQ'2+Q'iQ2-9^''Qi-Q2-ie''''''^''Qii,'Q2u') ■ (10) 

All vectors appearing in the above equation are 4-vectors, and T2 is obtained from Tl with the interchange of indices 
(5 A, 7 I/). The neutrino phase space integrals are easily performed by introducing a delta function as 

1 = J d^Q d\Q - Qi - Q2) , (11) 

and then using Lenard's identity 



N, 



J l^i;^Qi«Q2/3<54(Q-Qi-Q2) = ^(Q25„^ + 2Q„Q;3)e(Qo)e(Q2_q2) . (12) 

Before performing further contractions, we choose to retain only the transverse part of the photon propagator in 
eqn. (8). We will justify this a posteriori when the contribution from Landau damping is found to be much larger than 
naive screening [8]. We focus on this effect, which is expected to be important due to the exchange of soft photons. 
The remaining tensor contractions in eqn. (8) finally yield 

Q^^-Br = A J dQoQtP , (13) 
(27r)i2 n p%^ ' ^ ^ 
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nF{E^)nF{E^)hF{E^)fiF{E^)5\Vi - V,)6{Qo + E, - Ei) 

{|P2 xk|2 (13£;iE3-Pi-P3)} • (15) 



where boldface quantities denote 3-vcctors. The momentum integrals over d'^p2,d'^p^ arc performed first using the 
techniques described in [10]. Since the electron mass rUe pf^, the electron Fermi momentum, the electron energies 
Ei may be traded for the momentum magnitudes pi at the Fermi surface. Then P in eqn. (15) becomes 
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Next, the dpi,dp4 integrals are performed, which fixes the hmits on the dk and duj integrations. At this stage, it is 
convenient to perform the d^k integral, noting that sin^^o is a function of to, k. Explicitly, this integral is of the form 



r^'^' sin^^o 



ih= [ dk ^ , 7;\;", , (18) 



(^2 _ fc2)2 + 



M 



where Zi ^ T is the lower integration limit as determined by the G-functions involving k. The evaluation of this 
integral is detailed in appendix A. The result is 



--tan"^' ^ 



(19) 



Wc restrict our temperature and density range such that this approximation makes only a 10% error in the final 
analytic expression for the emissivity, which is acceptable. Using this approximation, P reduces to 

P = l§l!^«T^7V(a), (20) 
N{a)= Tdnn^e-^'m); = g , a=^ . (21) 



where I{^) is a smooth function of Cl obtained upon integration over uj (sec Appendix B for details). Substituting the 
resultant expression for P from eqn.(20) in eqn.(13), we are left with an overall integral over Qq, which is numerically 
tractable. The final result for the neutrino emissivity reads (we have chosen natural units ?i = c = 1; a factor of h^^c^ 
should be inserted in the denominator to recover "standard" units) 



Q,-.,= (^l ''"'df^'-i ma) . (22, 



""D 

Substituting the values of the physical constants, and defining Tg as the temperature in units of lO^K, nio as the 
electron chemical potential in units of 10 MeV, we obtain 

Qee-Br = 4.6 X lO^*^ fijo N{a) erg cm'^ g-^ . (23) 

where we also introduced a normahzed N{a) = N{a) /212 such that iV(Q:) < 1 (See Figure 2). The emissivity is plotted 
as a function of temperature for an electron chemical potential fig = 30 MeV in Fig. 3. For comparison, the emissivity 
as computed from eqn.(51) of ref. [8], which assumed Debye screening of both the longitudinal and transverse parts 
of the photon propagator, is also shown. The inclusion of Landau damping has modified the temperature; and density 
dependence of the emissivity from the case of pure screening. In our case, the emissivity scales as /Xg while in the 
screened case, it behaves as He T^. Note that in our result, the quantity A''(a) is also temperature dependent, so that 
a slight curvature in the log-plot of emissivity versus temperature is to be expected. In case of [8], a similar numerical 
integral appears but is temperature-independent and their result yields a straight line. The Tg dependence implies 
a slower fall-off of the emissivity with decreasing temperature, reflecting the increased importance of soft photon 
exchange at lower temperatures. This is due to the fact that electrons have smaller thermal energies, and cannot 
scatter far from the Fermi surface. The soft exchanges are not screened, and this leads to the enhancement of the 
emissivity, when compared to previous results which did not incorporate the Landau damping. At high temperatures 
(Tg » 10 and beyond), the electrons start to become non-degenerate, and our approximations would break down. 
Nevertheless, the convergence of our result to that of ref. [8] with increasing temperature is expected since both must 
approach the emissivity computed in the non-degenerate case [16]. 
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FIG. 2. The normalized N{a) = N{a)/212 from eqn.(21) 
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Emissivity vs Temperature 

()A =30 MeV) 



' " 1 1 ' 




1 




— Our result 






-- Ref. [8] 


y 


^ — 

- ^ 

— y 
y 


^^^^ ^ 

y 

y 

y 

y 


y 

y — 

y 




10'^ 10"' 


lo" 


10' 



Tg = (T/10 K) 

FIG. 3. The neutrino emissivity (on a log scale) from eqn.(23) of this paper and eqn.(51) of [8]. The chemical potential is 
fixed at p,^ = 30 MeV. 



There are other sources of neutrino emission from a QED plasma, such as the plasmon process [17,18], where an 
in-medium photon dressed by the electron plasma can decay into vv pairs, and the plasma photo-neutrino process [19] 

-I-7 ^ -|- vv. Formally, the plasmon process can interfere with vv bremsstrahlung when the exchanged photon in 
electron-electron collisions is almost on-shell. This is because the lifetime of such a photon is very long and it can be 
regarded as a valid in-medium excitation just like a plasmon. However, these two contributions can be disentangled 
when the temperature and density are such that only one of them dominates [20] . Studies of the plasmon process 
have shown that for a degenerate electron plasma at low temperatures such that T <C ti^p, where (jjp is the plasma 
frequency, the total emissivity acquires an exponential factor e~"^/"^. For example, the emissivity from the decay of 
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the transverse mode (photon) is ^ [18] 

^ ^^^ Jl^]^^/^T^/^e--^/T . (24) 
The plasma frequency for a relativistic degenerate plasma is a function of the electron chemical potential 

For u)p ~ 0{MeV) and for our temperature regime (T < 1 MeV), the plasmon decay rate and the plasma photo- 
neutrino rate are strongly suppressed, while the bremsstrahlung rate is not. In the following sections, we apply our 
main result, eqn.(23), to explore and compare neutrino rates in a variety of astrophysical environments. 

IV. COMPARISON TO NEUTRINO RATES IN STRANGE STARS 

Self-bound quark matter made of up, down and strange quarks can exist at high density, as evidenced by model 
calculations [21]. An entire star made of this phase is known as a strange star. Having obtained a modified in- 

mcdium neutrino radiation rate, wc would like to identify the physical conditions under which neutrino emission from 
relativistic electron-electron scattering, and quark-quark scattering, can affect strange star cooling. 

A. The strange star surface 

In case the star's surface is not covered by baryonic matter, it has a bare quark surface covered by a thin layer 
of degenerate electrons [22], termed the electrosphere. It has recently been suggested that the bare quark surface 
is not as sharp as previously thought but may rather consist of a crust of quark nuggets embedded in a degenerate 
electron gas [23]. These strange nuggets have larger electric charge than nuclei in a nuclear crust, but are very 
sparsely distributed, so that neutrino luminosity from electron-nugget bremsstrahlung will still be less than that from 
electron-electron bremsstrahlung, as calculated in this paper. Another intriguing scenario develops if the underlying 
quark matter is in a gapped phase, since the neutrino emissivity from quarks is severely attenuated at temperatures 
T <^ Tc where Tc is the critical temperature for pairing. The specific heat is similarly affected if electrons are not 
admixed with quark matter (as is the case with the color-flavor- locked or CFL phase [24]). Despite the enforced 
charge neutrality of the CFL phase in bulk, surface effects still lead to the formation of a degenerate electron plasma 
at the star's surface [25] and a possible nugget envelope along the lines of [23] may be present. Neutrinos from 
electron-electron bremsstrahlung in the envelope will then be the dominant energy sink for a strange star composed 
of stable CFL matter all the way up to its surface. In addition, the superfluid nature of the CFL phase implies a 
large thermal conductivity [26], so that neutrino cooling is determined by the neutrino emission properties of the 
surface [27] . Below, we examine the dominant sources for neutrino emission in the bulk of the strange star, since the 
true phase of quark matter at supra-nuclear density remains as yet unknown, and make a comparison to the surface 
or envelope emission rates. 

B. Ungapped quark matter 

If quarks of all flavors are ungapped, the direct urea process (weak equilibrium between d and u quarks) is dominant, 
with [28] 

Qdurca = 8.8 X 10^6 ae ( —) Y^'^'tI erg cm^^^-i ^ (26) 



^The emissivity can be decomposed into longitudinal and transverse parts for the vector channel, and an axial channel. They 
are all exponentially suppressed as exp(— wp/T) though the pre-factors are different [18]. 
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With typical values ac = 0.1, baryon density ub = 2no (no = 0.16 fm~^ is the saturation density of nuclear matter) 
and electron fraction = 10~^, a comparison of eqn.(26) with eqn.(23) shows that the former dominates by more 
than 8 orders of magnitude for the entire temperature range under consideration. With the inclusion of the non- Fermi 
liquid effects in ungapped quark matter, the emissivity is [14] 

2 

erg cm^^s"^ . (27) 

Other neutrino emission processes from ungapped quark matter, such as from quark-quark scattering, also exceed the 
rate computed here. Therefore, the presence of ungapped quark matter below the surface obviates a consideration of 
neutrino emission from the surface, at least as regards the cooling behaviour of the star. 

C. Gapped quark matter 

In the CFL phase, all 9 quarks are gapped, and the expectation is that neutrino emission processes involving free 
quarks are suppressed by the energy cost to be paid in exciting them above the gap A 3> T. This is not the complete 
story, however, since the CFL phase supports Goldstone modes (8 massive and 1 massless) that require much less 
energy to excite. Importantly, they can all couple to neutrinos, and contribute to the emissivity. Among them, the 
most important processes are the electroweak decay of the kaon, and of the exactly massless H boson through a 
momentum dependent coupling to neutrinos [29]. The latter has a very weak temperature dependence [27] and can 
be effectively ignored until very late times in the star's history. A comparison of the neutrino rate from kaon decay 
(eqn.(30) in ref. [27]) with eqn(23) is shown in Fig. 4, revealing that surface neutrino emission from the electrosphere 
starts to dominate for temperatures less than 8 x 10®K or so. The reason for the sharp fall-off of the neutrino emissivity 
in the CFL phase is the Boltzmann factor of exp(— mcs/T) where mcB, the mass of the Goldstone boson is typically 
10 MeV. 
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FIG. 4. Comparison of neutrino emissivity (on a log scale) from eqn.(23)(solid) with eqn.(30) of ref. [27] (dashed) , and 
eqn.(31) of ref. [36] (dash-dotted) . The quark chemical potential is fixed at Hq — 500 MeV and the electron chemical potential 
at fie = 30 MeV. 

Stresses placed on the CFL phase by virtue of a finite strange quark mass can trigger the formation of new ground 
states. For example, the lightest massive boson can Bose-condense, leading to neutral kaon condensation (the CFLK" 
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phase) [30-32]. While the decreasing mass can initiaUy increase the cmissivity from the decay ^, phase space 

considerations eventually shut off neutrino emission. The dominant rates come from processes such as —> e^Ve 
(and the analog for the K^) and K^e^ v^, which can occur due to the dispersion relation turning space-like beyond 
a certain momentum. The associated neutrino rates are enhanced over the CFL rates, but still suffer from exponential 
suppression due to Boltzmann factors, and cannot compete with the surface emission when the temperature drops 
below 0.1 MeV or so. 

The CFL and CFLK" phases arc; the; likely ground state of dense quark matter only at extremely high quark chemical 
potential. At densities relevant to neutron stars, with /Xg ~ 500 MeV or less, and with the physical requirements of 
charge and color neutrality, the pairing pattern can be quite complex involving phases with gapless modes for certain 
quark quasiparticlcs [33 35]. This opens up the possibility of a direct urea process that is not exponentially penalized 
even though the gap is non-zero, and the neutrino emissivity from gapless CFL matter (gCFL) was calculated in [36]. 
To facilitate a comparison of the rates, we choose (as in [36]), fiq = 500 MeV, m1/^,q = 50, and a gap A = 25 
MeV. From Fig. 4, it is clear that this is a powerful channel for neutrino emission at all temperatures relevant to the 
long-term coohng of the star. It is a sensitive function of ml/ ji, but the values used here are likely close to those in 
nature. Its effect on overall cooling only shows up in the photon- dominated cooling epoch however, due to the large 
specific heat of the gapless phase [36] . 

Several other color superconducting phases at intermediate densities relevant to neutron stars have been considered 
in the literature [37-39] which we do not address individually here. The question as to which is the preferred stable 
state of quark matter at intermediate densities and physical strange quark mass remains an open one at this time. 
We may mention here that a crystalline (LOFF) phase may exist at large values of the strange quark mass [40,41] and 
the neutrino emission from such a phase remains to be computed within the framework of the low-energy effective 
Lagrangian for this phase [42] . 

The enhancement of neutrino emission due to the inclusion of soft exchanges in electron-electron collisions can 
also be extended to quark-quark collisions, though the issue is much more subtle. Neutrino bremsstrahlung rates 
from quark-quark collisions have been previously computed in [28,43,44], although without damping in the transverse 
channel taken into account. A treatment of quark-quark scattering and transport properties in a degenerate quark- 
gluon plasma including this effect was undertaken in [10]. Having performed an explicit computation for electron- 
electron scattering, we may expect a simple extension to neutrino bremsstrahlung in quark-quark scatterings, and for 
different quark phases. For ungapped quark matter with chemical potential /z for a particular flavor, the emissivity 
from eqn.(22) applies to scattering between like-flavor quarks, with an extra factor of 4/3 from color-counting and 
averaging, and with the replacements e —> g for the coupling constant, /ig ^q, c^g Cq. The weak couplings 
for quarks Cq are listed in Table 1 of [45]. The Debye mass is m|, = acHq/n. However, as pointed out in section 
II, at energies (temperatures) u; ~ m e~^'^/^"* ~ 10s of keVs, non-Fermi liquid effects can render a perturbative 
analysis invalid, leading to logarithmic enhancements in the quark sclf-cncrgy and the neutrino rate in ungapped 
quark matter [13,14]. Already for Tg < 0.1, non-Fermi liquid effects are expected to set in. Since we integrate over all 
exchange energies, a more accurate evaluation of bremsstrahlung neutrino rates should include the non- Fermi liquid 
effects. In this respect, the neutrino rate computed here will be somewhat modified for quark matter, but only very 
slightly so for a QED plasma. Further complications arise if quark matter is superconducting. For gapped quark 
matter, such enhancements as discussed here need only be considered when an ungapped quark species exists, and 
interacts via a gluon whose Meissner mass (but not the Debye mass) vanishes. In the CFL phase, all quarks arc 
gapped, and all transverse gluons acquire a Meissner mass [46]. In the gapless phase, the issue of Meissner masses 
is far from clear [48-50]. However, a residual U{1) symmetry in this phase implies that for a condensate that is 
neutral under the associated conserved charge (denoted Q), the Meissner mass vanishes. Since some quark species 
carry Q charge however, the gCFL phase is not a Q insulator, and the Q photon (like the other gauge bosons in 
this phase) is likely screened. Finally, in the LOFF phase ^, whose appearance may be linked to the chromomagnetic 



^At first sight, sucli a process violates helicity conservation, but the presence of a preferred frame implies that the decay is 
indeed allowed in a frame moving with respect to the medium [27]. 

''The two-flavor superconductor, i.e, the 2SC phase admits ungapped quarks and unscreened gluons, but is energetically 
disfavored as a homogeneous phase in most studies [39,47]. 

® An alternate phase that can arise at intermediate density is the Overhauser phase, which like the LOFF phase, is anisotropic; 
however it seems not to be favored energetically unless the number of colors is large [51] or the number of space dimensions is 
reduced [52]. 
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instability of the gaplcss phases, the 5 gluons corresponding to the unbroken SU{2)c group acquire Debye masses 
due to the existence of blocking regions on the Fermi sphere, while the Meissner masses vanish [53]. Since u and d 
quarks of one color ( "blue" ) as well as strange quarks of all colors (presumably too heavy to participate in pairing) 
remain unpaired, unscreened gluon exchange among them can lead to enhanced neutrino bremsstrahlung emissivities, 
although the dominant neutrino emission process would be the modified urea among ungapped light quarks (eg. 
d + d^d + u + e~+Pe)^- 

V. COMPARISON TO NEUTRINO RATES IN BLACK-HOLE ACCRETION TORI 

Extreme conditions in which neutrino emission is crucial for the system evolution are also found in accretion tori 

around black-holes. Such tori are essential ingredients of most Gamma-Ray Burst models, cither in the collapsar 
scenario for long bursts [54] or in the neutron star-neutron star or neutron star-black hole merger scenarios for short 
bursts (see, e.g., [55] and [56]). With central densities above 10^^ g cm~^ and temperatures in the range of 1 to 10 MeV, 

nucleons form an almost ideal gas while electrons are mildly degenerate. Considering, for illustration, comparable 

1 /3 

fractions of neutrons and protons, one has for the electron chemical potential /x ~ 10 p^Q MeV. In the outer part of 
the torus, electrons are non-degenerate and thus positrons are present. 

A reliable comparison of the various neutrino emission processes can only be performed through numerical calcu- 
lations due to the wide range of physical conditions present in such tori and their rapid time variation, of the order 
of tens of milliseconds. We only present here some estimates based on the physical parameters that were employed 
in the figures of [56]. The dominant neutrino cooling process in these models of accretion tori is electron capture on 
protons [57], i.e., e~ + p ^ n + Vg, which in the degenerate electron limit, has a total emissivity 

Q^ll, = 4.2 X 10^1 M?o erg cm'^ s'^ , (28) 

while eqn.(23) gives 

Qec-Br = 4.6 X lO^" T^, fijo N (a) erg cm'^ s'^ , (29) 

showing that neutrino bremsstrahlung from e-e scattering can make a significant contribution in the "low" density 
and temperature regions of the torus, and even become the dominant emission mechanism. Maybe more important 
than energy loss is the production of antineutrinos since it is often considered that the jet of the GRB is triggered 
by pair production from u-V annihilation along the black-hole rotation axis. Neutrinos are copiously produced by 
electron capture but antineutrinos come essentially from the outer part of the torus, where T 10^° K and p < 10^° 
g cm~^, and are produced by bremmsstrahlung in nucleon-nucleon collisions, i.e., N + N' ^ N + N' + v + V, with a 
rate [58] 

QNN-Br = 4.7 X 10^1 Tlo^ pIo erg cm'^ s'^ , (30) 

where piQ is the density in lO'^'^ g/cc. Other important contributions arise from the outermost regions, at densities p <C 
IqW g (j]jj-3 -^yiigj-g electrons are non-degenerate and positrons are also present, from the process of pair annihilation, 
e~ + e'^ ^ u + V, with a rate 

Qpair = 4.8 X 10^4 erg cm'^ g-^ , (31) 

and positron capture on neutrons, n -|- e+ ^ p + 1'^, with a rate 

Qcap = 4.6 X 10^4 Tfo P7 erg cm'^ s'^ , (32) 

However, neutrinos from e-e scattering which, in contradistinction to processes involving positrons, come from within 
the whole torus volume, can be expected to make by far the dominant contribution to antineutrino production. 
Inclusion of this process in future numerical simulations will allow an assessment of its importance. 



^The quark direct urea process d ^ u + e + Ve would be suppressed due to the mismatch 5p = — fJ-u being large in this 
phase. 
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FIG. 5. Relevance of the e-e bremsstrahlung process in neutron stars. Left panel: comparison with the electron-ion 

bremsstrahlung neutrino emission in the crust, at three typical temperatures. Right panel: comparison with the modified 
Urea ("MUrca"), sum of the neutron and proton branches, and the proton-proton bremsstrahlung process in the core at the 
same temperature. 



VI. COMPARISON TO NEUTRINO RATES IN NEUTRON STARS 



Neutrino pair bremsstrahlung in electron-electron collisions can also be important in neutron stars and we discuss 
this in the present section. In the star's crust, the major neutrino emission processes arc pair bremsstrahlung processes 
for electron- ion scattering while in the neutron ^Sq superfluid regime, it is neutrino pair emission by breaking and 
formation of Cooper pairs at a temperature just below Tc, the critical temperature for pairing. The Cooper pair process 
is by far the dominant one but only acts efficiently in a narrow density region where the temperature T ^ 0.2 — 0.9 Tc, 
because of the density dependence of Tc, while the e-ion and e-e bremsstrahlung processes act in the whole crust. We 
plot in the left panel of Fig. 5 these last two processes: the emissivity of neutrinos in e-e collisions can dominate only 
at temperatures < 10^ K and the highest density regime where the "pasta" phase is present. One must emphasize that 
this "pasta" region actually comprises almost half of the mass of the crust and, moreover, neutrino bremsstrahlung in 
e-ion collisions is expected to be actually much weaker than shown in the figure [59] , by up to two orders of magnitude, 
due to the lower dimensionality of the nuclei (which are rods and plates instead of spheres). 

In the core of the neutron star, the modified Urea and nucleon bremsstrahlung processes are usually considered as 
the main neutrino emission processes in the absence of any enhanced neutrino emission process such as, e.g., the direct 
Urea. We plot in the right panel of Fig. 5 the emissivities of the modified Urea process as well as the bremsstrahlung 
rates in p-p (with rates taken from [60]) and e-e collisions: naturally, the modified Urea process is the dominant one. 
However, in regions where neutrons arc superfluid and T <^ Tc the modified Urea process is strongly suppressed, as is 
neutrino pair emission by the Cooper pair formation process and n-n and n-p associated bremsstrahlung: in such a case 
only the bremsstrahlung from p-p and e-e scatterings is effective and the latter one is dominant. Notice that in case 
of proton pairing, the bremsstrahlung rate from p-p scattering will become suppressed but that from e-e will also be 
significantly reduced because in a superconductor static magnetic fields are screened by the Meissner effect. Photons 
acquire a Meissner mass and electromagnetic interactions are thus fully screened instead of having their transverse 
part only Landau-damped. In such a case the calculation of [8] should be applicable for the bremsstrahlung emissivity 
from e-e scattering, since these authors assumed screening of the whole propagator by a Debye mass niD (instead of 
mjj and rriM)- Again, this process turns out to be the dominant neutrino emission process, although quite inefficient. 
In the outer part of the core, protons are certainly paired with a Tc of the order of 3 x 10^ K at most [61], while they 
are likely to be in a normal state at higher densities; the possibility that neutrons may pair at these densities is still 
very controversial and predictions range from possibly quite high values of Tc [62] to almost zero [63] . 



11 



VII. CONCLUSIONS AND OUTLOOK 



We have performed an explicit computation of the neutrino emissivity from neutrino bremsstrahlung in electron- 
electron collisions in a rclativistic degenerate plasma within an astrophysical context. This process has been studied 
before in the context of neutrino cooling from neutron star crusts, and in gapped superfluids. In this work, we have re- 
investigated this process in different astrophysical settings, with closer attention to the dynamics of the screened photon 
propagator that mediates electron-electron collisions. Compared to previous results [8], we find an enhancement when 
Landau damping in the transverse part of the photon propagator is taken into account. If we assume simple screening 
in both the longitudinal and transverse channels of photon exchange, we recover the results of [8]. Our principal 
results are eqn.(23) and the comparisons in Fig. 3, Fig. 4 and Fig. 5. The inclusion of Landau damping implies a /ig T"^ 
dependence, instead of the jip dependence as expected from screening arguments. This leads to an enhancement in 
the neutrino rate which increases with dc;creasing temperature, suggesting that e-e bremstrahlung is indeed an efficient 
cooling process in a certain window of temperature and density. This conclusion could be of particular importance 
to cooling of strange stars when the interior of the strange star is in a color superconducting phase, such as the CFL 
phase, which has only gapped quark excitations. The neutrino emission from the interior is then suppressed by several 
orders of magnitude as the temperature falls below the critical temperature for pairing (T^ ~ 50 MeV). The overall 
cooling of the star can be affected due to surface effects as shown in [64,65]. Inclusion of the enhanced neutrino rate 
in the computations of [64] would be a natural extension of this work, and taken in conjunction with the photon 
cooling from the strange star surface (which can also influence the light curve at early times) serve to present a more 
complete picture of strange star cooling. 

The quark matter underlying the strange star surface can exist in any one of many phases, and neutrino emission 
rates have been computed in only a few of them. Of particular interest are the gapless and crystalline phases, which 
will lead to temperature dependences that are modified from the usual forms in ungappcd rjuark matter. These phases 
also have unique dispersion relations for certain quark quasiparticles, and consequently, a specific heat per unit volume 
that is also different from ungapped quark matter. These two factors imply a change in the stellar cooling curve that 
can be confronted by observations. While we have focused on computing a particular process, our ultimate goal in 
this direction is to provide more accurate neutrino rates, including the possibility of various quark phases, in order 
to make accurate predictions for the cooling history of a strange star. This task is complicated by the appearance of 
non- Fermi liquid behaviour, and the breakdown of a perturbative treatment at temperatures relevant to the long-term 
cooling of quark matter in the stellar interior. 

Finally, we have established the relevance of neutrino bremsstrahlung in e-e collisions for black-hole accretion tori 
and crustal layers of ordinary neutron stars. The GRB jet may be triggered by vD annihilation, while its subsequent 
cooling also depends on the neutrino flux, and the e-e process is important for both. Neutron star cooling can also 
be affected by this process at low temperatures in the pasta-phase of the crust or when neutrons are superfluid in the 
core. Our results thus find wide-ranging applications for neutron stars, GRB's and strange stars. Increasing data on 
neutron star cooling, GRB afterglow signatures and theoretical advances on observable signatures of strange stars [66] 
all suggest that a renewed effort to pin down neutrino rates for astrophysics with improved accuracy is both timely 
and interesting. 
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APPENDIX A: JsT-INTEGRATION 



In eqn.(18), a change of variable to x = uj/k gives 




(A2) 



(Al) 



12 



For the physical conditions of interest, a <C 1 and 6 1, which enables us to use a much simpler form 

= ^"irT6^ • ^^^^ 

The integrand is not very different than in eqn.(Al) for the entire range of integration (note that Zi ^ T so that 
w/Zi ~ 1), as evident from Fig. 6 below. 



0.08 r 



0.06 



0.04 



0.02- 




FIG. 6. The integrands from eqn.(Al) and eqn.(A3) labeled g{x) and f{x) respectively. The temperature is Tg = 10 
(T = 10^°K) and jjno = 3 (/Xe = 30 MeV). 



The integral in eqn.(A3) is straightforward and yields 
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APPENDIX B: NUMERICAL EVALUATION OF N{a) 



The integral in eqn.(21) is 
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(B2) 



This integral is well-behaved for all finite values of fl, even as w — > 0. At O = 0, the integrand diverges at a) = 0, but 
the overall cmissivity is clearly zero. Thus, a numerical evaluation of N{a) from eqn.(21) with the above expression 
for /(fi) is feasible. The rapid exponential fall-off of the integrand with respect to u and Cl implies a rapid convergence 
and high degree of accuracy in the final result. The integral N{a) is evaluated by Gauss-Legendre quadratures and 
the normalized N that appears in eqn. (23), is plotted in Fig.2. 
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